[!
MrBayes settings for the best-fit model (GTR+I+G) selected by AIC in MrModeltest 2.2

BEGIN MRBAYES;
	Prset statefreqpr=dirichlet(1,1,1,1);
	Lset  nst=6  rates=invgamma;
END;
!
MrBayes settings for the best-fit model (HKY+I+G) selected by AIC in MrModeltest 2.2

BEGIN MRBAYES;
	Prset statefreqpr=dirichlet(1,1,1,1);
	Lset  nst=2  rates=invgamma;
END;
]

begin mrbayes;

outgroup russicum;

end;



begin mrbayes;		
       		Lset  applyto=(1) nst=6  rates=invgamma;
      		Prset applyto=(1) statefreqpr=fixed(equal);
		Lset  applyto=(2) nst=2  rates=invgamma;
		Prset applyto=(2) statefreqpr=dirichlet(1,1,1,1);
     	        unlink shape = (all);	
		unlink tratio = (all);	
		unlink revmat = (all);	
		unlink statefreq = (all);
		unlink pinvar = (all);
		unlink omega = (all);
	    set autoclose=yes;	
		mcmcp ngen=150000 printfreq=1000 samplefreq=100 diagnfreq=1000 nchains=4 savebrlens=yes  filename=BIGRUN1;	
		mcmc;	
		sumt filename=BIGRUN1  burnin=100 conformat=simple contype=halfcompat;
end;